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Abstract 



This work examines the use of two-way training in multiple-input multiple-output (MIMO) 
wireless systems to discriminate the channel estimation performances between a legitimate 
receiver (LR) and an unauthorized receiver (UR). This thesis extends upon the previously 
proposed discriminatory channel estimation (DCE) scheme that allows only the transmitter 
to send training signals. The goal of DCE is to minimize the channel estimation error at 
LR while requiring the channel estimation error at UR to remain beyond a certain level. If 
the training signal i= sent only by tl>e ttaiSl^ performance discrimination bet«en 
LR and UR will be limited since the training signals help both receivers perform estimates 
of their downlink channels. In this work, we consider, instead the two-way training method- 
ology that allows both the transmitter and LR to send training signals. In this case, the 
training signal sent by LR helps the transmitter obtain knowledge of the transmitter-to-LR 
channel, but does not help UR estimate its downlink channel (i.e., the transmitter-to-UR 
channel). With transmitter knowledge of the estimated transmitter-to-LR channel, artificial 
noise (AN) can then be embedded in the null space of the transmitter-to-LR channel to 
disrupt UR's channel estimation without severely degrading the channel estimation at LR. 
Based on these ideas, two-way DCE training schemes are developed for both reciprocal and 
non-reciprocal channels. The optimal power allocation between training and AN signals is 
devised under both average and individual power constraints. Numerical results are provided 
to demonstrate the efficacy of the proposed two-way DCE training schemes. 

i 



Contents 



Abstract 



Contents 



1 Introduction 




3 Two- Way Training Strategy for. Reciprocal Channels 




4 Optimal Power Allocation for DCE in Reciprocal Channels 

4.1 Channel Estimation Performance at Transmitter 

4.2 Channel Estimation Performance at LR and UR 

4.3 Optimal Power Allocation between Training and AN Signals . . . 

5 Two- Way Training Strategy for Non-reciprocal Channels 

6 Optimal Power Allocation for DCE in Nonreciprocal Channels 

6.1 Channel Estimation Performance at Transmitter 

6.2 Channel Estimation Performance at LR and UR 



6.3 Optimal Power Allocation between Training and AN Signals 

Numerical Results and Discussions 

Conclusion 

Appendix 

9.1 Proof of Proposition I 

9.2 Monomial Approximation and Condensation Method for the Problem in fl6.33p 




iii 



List of Figures 



2.1 Diagram of a wireless MIMO system consisting of a transmitter, a legitimate 
receiver (LR) and an unauthorized receiver (UR) 



7.1 NMSE performance versus forward training Tp for Pt = 30 dB and Pl = 20 dB. |26 

7.2 Power allocation among reverse and forward training powers Er/tr-, Sf/tf 
and AN power (A^^ - NL)al 

7.3 Power allocation among the training powers So/Nt, £i/Nt £2/^1 Ss/Nt and 



27 



AN power (A^ - A^)(t^. . \k 

7.4 NMSE performance of the proposed DCE scheme for the reciprocal and non- 
reciprocal case 

7.5 SER performance of the LR and UR in an OSTBC system with the channel 
estimates obtained by the proposed DCE scheme 



28 



30 



31 



iv 



Chapter 1 
Introduction 



Secrecy in wireless communications has been an important problem over the years due to the 
broadcast nature of the wireless medium. In the past, these issues have mostly been addressed 
using cryptography in the application layer. However, recent studies on information-theoretic 
secrecy provide an alternative to achieve these tasks through coding and modulation in the 
physical layer. In the context of physical layer secrecy, one is often interested in deriving 
the so-called secrecy capacity, which is the rate achievable with vanishing error probability 
at the legitimate receiver (LR) and vanishing equivocation rate at the unauthorized receiver 
(UR). In particular, secrecy capacity has been derived for single- input single-output (SISO) 
systems in [T] and for multiple-input multiple-output (MIMO) systems in |2j. The results 
show that secrecy capacity can generally be increased by enlarging the difference between the 
effective channel qualities of LR and UR. While most works on physical layer secrecy focus 
on optimal coding schemes to achieve secrecy capacity in the data transmission phase, our 
goal is to exploit signal processing methods to enlarge the differences between the quality of 
the two channels. In particular, this is done from a channel estimation aspect, following the 
so-called discriminatory channel estimation (DCE) methodology proposed previously in [3]. 

Specifically, DCE is a training strategy that utilizes artificial noise (AN) to disrupt UR's 
reception while sending training signals to LR. In this case, AN must be placed in the null 
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space of the transmitter-to-LR channel to minimize its effect on LR. However, this requires 
transmitter knowledge of the channel, which is typically obtained through feedback from 
LR. In the original DCE scheme, only the transmitter is allowed to send training signals. 
In this case, increasing the training power helps improve the channel estimate at LR and 
allows for more effective use of AN at the transmitter. However, this also helps UR obtain 
a better channel estimate and, thus, the amount of power used for training must be limited. 
To improve the channel estimate at LR while confining the performance of UR to a certain 
level, multiple stages of feedback and retraining must be employed. In this case, training 
and AN power is increased as the transmitter knowledge of the channel improves through 
multiple stages. Yet, the training overhead and complexity required to optimize training 
over multiple stages limits its application in practice. 

The main contribution of this thesis is to propose new and efficient DCE schemes using 
the two-way training methodology. Here, 'tiraliSffig^^ignals will also be sent by LR and trans- 
mitter knowledge of the channel will be obtained by performing channel estimation at the 
transmitter. Notice that the original DCE scheme assumes that channel feedback with infi- 
nite resolution is provided from LR, which is not achievable in practice. When the channel is 
reciprocal, e.g., in time-division multiplexing (TDD) systems, the channel state information 
(CSI) can be obtained at the transmitter by sending pilot signals from the receiver. Two- 
way training schemes have been studied for conventional point-to-point links in [U [5l |6] to 
obtain the CSI at both the receiver and the transmitter without the use of feedback. In this 
work, we adopt the concept of two-way training into the design to increase the efficiency of 
the DCE scheme. In reciprocal channels, the proposed two-way DCE scheme uses reverse 
training to provide CSI at the transmitter and forward training with AN to achieve different 
channel estimation performances at LR and UR. When the channel is non-reciprocal, e.g., in 
frequency- division multiplexing (FDD) systems, the downlink and uplink channels between 
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the transmitter and LR would not be identical. In this case, an additional training phase 
is needed, where the transmitter first broadcasts a randomly generated signal to LR, which 
then echoes the signal back to the transmitter. The echoed signal contains information of 
both the downlink and uplink channels and can be combined with the reverse training signal 
to estimate the desired transmitter-to-LR channel {i.e., the downlink channel). Compared to 
the multi-stage feedback-and-retraining DCE scheme in [3], the proposed two-way training 
scheme drastically decreases the overall training overhead and design complexity. 

To optimize the performance of the proposed two-way training scheme, we derive the 
optimal power allocation between the training and AN by solving an optimization problem 
that aims to minimize the channel estimation error at the LR subject to a lower limit 
constraint on the channel estimation error at the UR. In the reciprocal case, the analytical 
result shows that the problem of finding the optimal training and AN powers reduces to 
a one- variable problem which can be solved by'^^t" simple line search. However, in the non- 
reciprocal case, the power allocation problem is not easily solved since the estimation error 
expression is much more complex. Therefore, we instead resort to an approximate solution 
by using the monomial approximation and Condensation method pi] in the field of geometric 
programming (GP). Numerical results show that the proposed DCE design can effectively 
discriminate the channel estimation and the data detection performances at the LR and UR. 

The remainder of the thesis is organized as follows. In Chapter [21 we first introduce the 
wireless MIMO system model considered in this work and provide a general description of 
the DCE scheme. For the case with reciprocal channels, the training strategy is described in 
Chapter E] and the optimal power allocation is derived in Chapter HJ Similarly, for the case 
with nonreciprocal channels, the training strategy is described in Chapter [5] and the optimal 
power allocation is derived in O Numerical results are provided in Chapter [7] and, finally, a 
conclusion is given in Chapter |8l 
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Chapter 2 
System Model 



Consider a wireless MIMO system that consists of a transmitter, a legitimate receiver (LR), 
and an unauthorized receiver (UR), as shown in Fig. 12.11 We assume that the transmitter, 
LR, and UR are equipped with Nt, and Nu antennas, respectively. The channels of LR 
and UR remain constant during one transmission block, which consists of a training phase 
and a data transmission phase. The goal is to prevent the UR to extract information from 
its received signal. Instead of focusing the data transmission, we propose to achieve this 
task from a channel estimation perspecti'Ve''to«f a&vise two-way training schemes following 
the DCE methodology that enables LR to perform an accurate estimate of the channel while 
disrupting the channel estimation performance at UR. 

Let the downlink channel from the transmitter to LR be denoted by G C^*^^^ and 
the uplink channel from LR to the transmitter be denoted by H^j G C^^^^*. In the following 
chapters, we consider separately two different channel models, i.e., the reciprocal channel 
model and the non-reciprocal channel model. In both cases, the proposed two-way training 
scheme for DCE can be generally divided into two steps as described below. 

Step I: The aim of Step I is to allow the transmitter to obtain an estimate of the 
downlink channel. Different from [3], where a noiseless feedback channel is required, we 
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Figure 2.1: Diagram of a wireless MIMO system consisting of a transmitter, a legitimate 
receiver (LR) and an unauthorized receiver (UR). 

allow the transmitter to estimate the downlink channel itself through the exchange of training 
signals between the transmitter and UR. This channel knowledge will be used in designing 
the forward training signal in order to discriminate the channel estimation performances 
between LR and UR. Different training strategies are.i'equired under different channel models 
to achieve the downlink channel estimsttion at the transmitter. Detailed descriptions for the 
reciprocal case and non- reciprocal case are given m^'Ghapters [3] and O respectively. 

Step II: After obtaining the downlink channel estimate in Step I, the transmitter next 
sends a training signal along with AN to degrade the channel estimation performance at UR. 
Specifically, by assuming that Nt > Nl, the forward training signal is given by 

where Ct G C^^^^' is the pilot matrix with Tr(Cf^Ct) = Nt, Vf is the forward training 
power, and Tp training length. For ease of notation, we define Sp — VpTp as the forward 
training energy. A G C^F^^^t-^^) jg ^^j-^g y^jsj matrix of which each entry is i.i.d. Gaussian 
with zero mean and variance o"^ and G C.'^t^(^t-^L) jg matrix whose column vectors 
form an orthonormal basis for the left null space of H^, that is, = Qin^_Nl)-xNl 
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the (Nt - Nl) by Nl zero matrix) and Nf Njj = lAr,_iv^. Notice from ([21]) that AN is 
superimposed on the training signal and placed in the left null space of to minimize its 
interference on LR. The received signals of the LR and UR are respectively given by 



where G G C'^t^^u jg ^\^q channel matrix from the transmitter to UR, and W G C^^^^^ and 
V G C'^^^^'^ are the additive white noise matrices at LR and UR, respectively. Each entry 
of G is assumed to be i.i.d. distributed with zero mean and variance equal to cr^. Elements 
of both W and V are assumed to be i.i.d. random variables with zero mean and variances 
respectively equal to and cr^. 

In the following chapters, we describe the training strategies and examine the correspond- 
ing channel estimation performances at the transmitter and the receivers during each stage 
of the process. The optimal power allocation betweeii training and AN signals are derived 
to achieve discrimination between the channel estimation performances at LR and UR. 





(2.3) 
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Chapter 3 



Two- Way Training Strategy for 
Reciprocal Channels 



In this chapter, we consider the case where the channel between the transmitter and the 
LR is reciprocal, which means that the downlink and uplink channels are symmetric. In 
this case, we define the downlink channel matrix as = H and the uplink channel matrix 
as Hu = H"^. With channel reciprocity, the transmitter can obtain an estimate of the 
downlink channel by taking the transpose of the channel matrix obtained through reverse 
training, i.e., training based on signals sent from^-bR to the transmitter. In this case, DCE is 
effectively achieved using only two stages, i.e., the reverse and the forward training phases. 
The operations are detailed below. 

Reverse Training: In the reverse training stage, LR first sends a training signal, denoted 
by X/^ G C^^^^^, to enable channel estimation at the transmitter. Specifically, the reverse 
training signal is given by 



where the pilot matrix Cl satisfies C^Cl = '^n^ (the Nl by A^^^ identity matrix), and Vr 
and tr represent the transmission power and training interval, respectively. For ease of use 
later, we define the reverse training as = VrTr. The received signal at the transmitter is 




(3.1) 
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given by 



Yt = X^H^ + W, 



(3.2) 



where each element of H is assumed to be independent and identically distributed (i.i.d.) 
random variable with zero mean and variance equal to aj^, and W G C^^^^' is the additive 
white noise matrix with each element having zero mean and variance a^. By the help of 
reverse training, the channel estimate of H, denoted by H, can be obtained at the transmitter 

Forward Training: In the forward training stage, the transmitter superimposes AN 
on top of the training signal to degrade the channel estimation performance at UR. With 
knowledge of the estimated downlink channel, i.e., H, AN can be placed in the left null space 
of H to minimize its interference on the LR. The forward training signal is given in (12.11) 
where the pilot matrix satisfies C^Ct = Iat^. And the received signals of LR and UR are 
respectively given in (12. 2p and (12. 3p . Note that the notation in (12. ip and (12.21) is replaced 



In the next chapter, we analyze the channel estimation performances of transmitter, LR 



criterion for channel estimation [7]. We then propose to judiciously allocate the training 
powers and the AN power in reverse and forward training, aiming at discriminating between 
the channel estimation performances of the LR and the UR. 




and UR by assuming that all of them employ linesEr minimum mean square error (LMMSE) 
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Chapter 4 



Optimal Power Allocation for DCE in 
Reciprocal Channels 



4.1 Channel Estimation Performance at Transmitter 

Due to channel reciprocity, the reverse training signals sent by LR allow the transmitter to 
obtain an estimate of the downlink channel by taking the transpose of its estimate of the 
uplink channel. By employing the LMMSE estimator, the estimate of the channel matrix H 
can be written as 

H = Kxf(4X,Xf + a^I.J-iY,f 

= H + AH (4.1) 

where AH G C^*^'^^ stands for the estimation error matrix. The covariance matrix of AH 
can be shown to be [7] 

E{AH(AH)''} = 7V.(^ + ^)"l„.. (4.2) 
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4.2 Channel Estimation Performance at LR and UR 



To analyze the channel estimation performance of LR, let us write (I2.2p as 

Yi = X —CtU - ANf AH + W = CH + W, 



(4.3) 



where C ^ J^Ct, W ^ 



AH + W and the first equality is due to N^H = 0. 



Denote the channel estimate at LR by H^,. The normalized mean squared error (NMSE) of 
Hi under LMMSE criterion can be shown to be [7] 

TrfE{(H-Hi)(H-Hi)^} 



NMSE, ^ 



Tr((Rj^^ + C^R^^C)"^ 



(4.4) 



where Tr(-) denotes the trace of a matrix, = Niajjl^^ and R,^ = E{WW^} is the 
covariance matrix of W. According to the .independence between A and W, the fact of 
and (14. 2p . it can be shown that 

R^= (E{||N|AH||2K + iVz.ay I, 

= Nl 

Substituting fOj) into f Oj) yields 



•-Tp 



2 I 2 

+ c^,„ 



I. 



(4.5) 



tr 



N, 



NtNL 



NMSE/ 



(Nt-NL) 



— + 



NtNL 
Sp/Nt 



(iv.-iv.)(4 



-1 



(4.6) 



The NMSE performance of the UR can be analyzed in a similar way. Specifically, one can 
show that the NMSE of estimating G at the UR is given by 



NMSEr; = {^ + 



Sp/Nt 



al {Nt-NL)alal + al 



(4.7) 
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4.3 Optimal Power Allocation between Training and 
AN Signals 



Observing from (14 .Gp and ( \A.7\i . the added AN in forward training can affect both the LR and 
UR's channel estimation performances. To optimize LR's channel estimation performance 
while preventing the UR from obtaining an accurate estimate of G, we propose to jointly 
optimize the reverse training energy S^, the forward training energy Sp and AN power 
by considering the following power allocation problem 



where we aim to minimize the LR's NMSE subject to a preset lower limit 7 on the UR's 
NMSE, under an average power constraint- Pa^e-r- l'^ote that the LR and the transmitter also 
have their own peak power constraints, i.e., Pl and Pt- 

Remark: In the DCE scheme, it is desirable to keep the forward training length as 
small as possible, i.e., equal to the number of transmit antennas. Observing from f l4.6p . 
( 14. 7p and the problem (14. 8 p and assuming the average energy constraint and the individual 
energy constraint on the transmitter and the LR are all fixed, as the forward training length 
increases, it needs more AN energy to meet the same lower limit value and thus less energy 
can be allocated to the training signal. Different from the receiver's noise of which the energy 
can be freely accumulated over time, it takes the system's resource to maintain the AN's 
power. 



min NMSEl 
s.t. NMSEc; > 7, 



(4.8) 



£r + £r + {Nt - NL)alTF < Pave{rR + Tf) 
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To make all constraints effective, we shall focus on the interesting case where 

max{PLTn, PtTp} < PaU^R + tf) < PlTr + PtTp. (4.9) 

Note that for the case where PaveijB. + ^f) > Pl^r + Pt'Tpi the average power constraint be- 
comes redundant and hence, the transmitter and the LR simply transmit with its maximum 
power. When PaveijR + tf) < PlTr and/or PaveijR + tf) < PtTp, one or both individual 
power constraints become redundant. The solution for this case can be easily obtained by 
following the derivations for the case of (14. 9 p 3. 

On the other hand, it should be noted that the preset value 7 should satisfy [3] 

1 mm{ PtTF, PaveiTR + '^f)}\ ^ 9 n 

^ ---^ R ni\ <7<4, (4.10) 

since the left-hand- side term is the minimum achievable NMSE of UR (when the transmitter 
does not use AN, i.e., al = 0), and the right-hand-side term stands for the worst NMSE 
performance of UR, respectively. For ease of tafter use, let us define 

7 = /-.-^,V.^?.>0. (4.11) 



Then the condition in (14.101) reduces to' ^^^'■'''■.^ 



< 7 < PtTF. (4.12) 

The power allocation problem in (14. 8p is a nonconvex optimization problem involving 
three variables {Sr,Sf,o''^)- However, it actually can be solved very efficiently. We show in 
Appendix 9.1 the following proposition for problem (14. 8p : 

Proposition 1. Consider the power allocation problem in (14. 8 p with both (14. 9 p and (I4.12p 
satisfied. If 



2 2 2 

f^ = NL{ ^ - ^ ) > min{P^r^, P^^rR + Tf) - 7} 



^The proposition to be given for the case of (|4.9I) also describes the sohition for the case of Pave {tr + tf) < 
Pltr and/or PaveiTR+Tp) < PfTp, by changing the condition in (|4.12p to < 7 < mm{PtTp, Pave{TR+Tp)} 
and setting the redundant individual power constraint (s) to infinity. 
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then the optimal {£ £ p ■, (^t) of f l4.8p is given by En* = 0, Sp = and (cr^)* = (i.e., 
no need of reverse training and no need of AN in forward training). On the other hand, if 
/i < mm{PLTR, PaveiTR + Tp) — l} , the optimal Sr of (14.81) can be obtained by solving the 
following one- dimensional problem 

^ w 

S.t max{0, /i, Pavei-TR + Tp) - PtTp} <£r< 
m\n{PLTpi, Pave{rR + Tp) - 7}, 



where 



and 



Tp + a^Gl al 



£p{£R) = l[^-a{8R) + iy (4.15) 
The optimal £p and a1 are given by £p = £p{£^) and (cTq)* = (jVt-jvi) • 

Proposition 1 implies that the solutions of problem (14. 8 p can be efficiently obtained by 
simple line search over a finite interval, when the condition in (I4.13P is fulfilled; otherwise, 
one can obtain a simple closed-form solution of £}^ = 0, £p = ^ and (cr^)* = 0. 
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Chapter 5 

Two- Way Training Strategy for 
Non-reciprocal Channels 

In this chapter, we consider the case of non-reciprocal channels, where the downlink and 
uplink channel matrices are asymmetric. In this case, the downlink channel cannot be 
directly inferred from the uplink channel;;, jireref ore, an additional training stage using an 
echoed signal (from transmitter to LR and back to the transmitter) is needed in order to 
obtain an estimate of the combined downlink and uplink channel. This additional stage 
is referred to as round-trip training. The; properSed two-way training method for DCE in 
non-reciprocal case is detailed below. 

Round-trip Training: In round-trip training, the transmitter first broadcasts a random 
signal then the LR will echo its received signal back. By the round-trip procedure, the echoed 
signal obtained at the transmitter contains a combined term of uplink channel and downlink 
channel. Then with the help of the following reverse training, the transmitter can obtain the 
downlink channel estimate. Specifically, the random signal sent by the transmitter is given 

by 

X,o = y^C,o, (5.1) 
where C^o G C^°^^' is the pilot matrix satisfying Tr(C^Cjo) = Nt, and Vq and Tq represent 
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the training power and training length, respectively. For ease of use later, we define the 
round-trip training energy as = VqTq. The received signal at the LR is given by 

Yio = X,oH, + Wo, (5.2) 

where each element of H^^ is assumed to be i.i.d. complex Gaussian random variable with 
zero mean and variance equal to aj^^ and Wq G C^"^^^ is the additive white Gaussian noise 
(AWGN) matrix with each entry having zero mean and variance cr^. Then the LR amplifies 
and forwards its received signal back to the transmitter. The echoed signal at the transmitter 
is given by 

Ya = aYMH„ + Wi (5.3) 
= aXioHdH„ + aWoH„ + Wi 

where each element of H„ is assumed to be i.i.d. complex Gaussian random variables with 
zero mean and variance crj^^, Wi G C^oxNt jglJig AWGN matrix at the transmitter with the 
power of each entry equal to a^. The Amplifying gain at the LR is given by 

(y = . - — — 

where Vi is the transmission power and £i = ViTq is the energy of the transmitted symbol. 
Since the random signal X^q is available at the transmitter, it is able to obtain the downlink 
channel estimate with a given uplink channel. We will see how a reverse training helps the 
transmitter to extract the knowledge of downlink channel H^. Note that the random signal 
X^o is unknown to both LR and UR, therefore the UR can not benefit from the round-trip 
training. 

Reverse Training: In reverse training, the LR sends a training signal G C^^^^^ 
to enable the uplink channel estimation at the transmitter. Specifically, the reverse training 
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signal is given by 



^L2 = (5.5) 



where Cl2 is the pilot matrix which satisfies Tt{C^2^L2) = Nl, and V2 and T2 is the 
transmission power and training interval of the LR. For simplicity, We define the reverse 
training energy as £2 = 'P2T2. The received signal at the transmitter is given by 

Y,2 = Xi2H„ + W2 (5.6) 

where W2 G C^^xNt jg ^j^g additive white noise matrix with each entry having zero mean and 
variance o"^. As the uplink channel estimate is given, the downlink channel estimate can be 
acquired from the echoed signal. 

Forward Training: In forward training, the transmitter sends AN along with the train- 
ing signal to discriminate the channel estimation performances between LR and UR. The 
specific description is stated in the Step II of Chapter |2J Note that we replace the subscript 
by 3 for notation consistency in this chapter, therefore the received signals at the LR and 
UR are replaced by *ts^-^a ui*^^-'-^x'' 

^3 TT , A tvtH 




Yl3 = \ CisH, + AN| H, + W3 (5.7) 
Yu3 = \^Ct,G + AN| G + V3 (5.8) 

and the forward training length is to substitute Tp- 

Due to the complicated nature of the two-way training, finding the optimal pilot struc- 
tures may be a difficult task, which could also be different for different objective functions, 
e.g. channel estimation error, bit error rate or ergodic capacity, etc. The practical intuition 
in choosing the pilot structure is 1) to reduce the channel estimation error and 2) to reduce 
the transmission overhead. In conventional channel estimation, the orthogonal structure 
was usually found to be good. Note that it may not be the optimal choice for the system 
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we are considering. By utilizing the orthogonal training signal, the performance of channel 
estimation is now determined by the training energy of each phase and one can keep the 
training length minimum if the training energy can be designed to reduce the channel esti- 
mation error. Besides, it is preferred to keep the training length small in the secrecy channel 
estimation according to the remark stated in Chapter HI Hence, in this work, we choose 
the minimum training length to be the number of transmit antenna, i.e., tq = ra = Nt and 
T2 = Nl. And we assume the unitary pilot data are used, that is C^Cfo = CtoC^ = '^Nt^ 
C^2^L2 = Cl2C^2 = and C^C^s = CtaC^ = Iat^. 
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Chapter 6 

Optimal Power Allocation for DCE in 
Nonreciprocal Channels 

In this chapter, we show how the transmitter can compute the downhnk channel estimate 
from the training signals and analyze the channel estimation performance at both LR and 
UR. We assume that the transmitter, LR and^-jSJtaliemploy the linear minimum mean square 
error (LMMSE) criterion for channel estimation [7j. Then, we examine the optimal power 
allocation between training and AN signals in this case and propose an efficient solution for 
this problem. " . ■ 

6.1 Channel Estimation Performance at Transmitter 

In this section , we show how the transmitter computes the downlink channel estimate from 
the reverse and round-trip training signals. Specifically, with the help of reverse training 
and by employing the LMMSE estimator, the estimate of the uplink channel can first be 
computed as [7] 

H„ = 4„Xf2(4„Xi2Xf2 + allN,r'^t2 = H„ + AH. (6.1) 
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where AH^ G C^^^^* is the estimation error matrix with correlation matrix given by 



E{(AHj^AH4 = iV^(— + -^1 I^,. (6.2 



1 ^2 



w 



With the uphnk channel estimate H„ being available at the transmitter, we can rewrite the 
echoed signal fl5.3p as 

Yn = aXioH,(H, - AH„) + aWo(H„ - AH„) + Wi (6.3) 
= aXioH,H„ + (-aWoH. - aXioHdAH, + aWoAH„ + Wi). 

To employ the LMMSE criterion for the downlink channel estimation at the transmitter, it 
is easier to empress (16. Sp in the vector form as 

yti = a(H^ ® Xio)hrf - a{AHj (^Xto)K + "(H^ ® ItvJwq - a(AH/ ® IjvJwq + wi (6.4) 

where the fact that vec{ABC) = {C^ ® Ayf&si^S) is used, y^i = vec(Y(i) is formed by 
stacking the columns of Y^i and so do h^^ = vec(Hrf), wq = vec(Wo), and wi = vec(Wi). 
As H„ is given at the transmitter, by the fact that H„ and AH„ are uncorrelated due to the 
orthogonality principle [7], the premise of CtoC^- — C^Cfo = Iatj and (16.20 . the LMMSE 
estimate of downlink channel h^^ and thus its matrix form are respectively given by 



+ (H:((HrH:)+/3u)"^x-)y, (6.5) 

^ h, + AVi (6.6) 

Hd, = A f ^ + ^fo^ti ((Hf H«) + P^n) Hf (6.7) 

^ Hrf + AH,,i (6.8) 



where 
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and Ahdt e C^t^^x^ is the estimation error vector at the transmitter. The correlation 
matrix of Ah^ ^ conditioned on a given H„ is given by 



E{Ah,,,(Ahrf,i)^|Hj 



Nt 



(6.10) 



Note that for differentiating from the downhnk channel estimate of the LR, we denote the 
downlink channel estimate of the transmitter as lid,t- The matrix consisting of the basis of 
left null space of Hd^t is replaced as ^ ■ 

6.2 Channel Estimation Performance at LR and UR 

In this section we analyze the channel estimation performance of the LR and UR. We first 
consider the channel estimation at the LR. Due to the fact that Hrf+ = the received 
signal of LR fl5.7p can be written as ,■■ //''-^f.^^^:j^i\- 



Yi3 = Q;,B|^.AN| AH,,, + W3. (6.11) 



where Ct3 = y j^Cts- To apply the LMMSE criterion for the downlink channel estimation 
of the LR, let us vectorize f l6.1ip as 

yL3 = {iNr. ® C«) hd - (liv^ ® AN^JAhd^t + W3 (6.12) 

where h^ ^ = vec(H(i) and W3 = vec(W3). Then the channel estimate of is given by 

hrf = Ch,yL,C-l^y^.jL3 (6.13) 

where 

C,,,,3 = E{h,yf3} = 4^ (I;v, ® Cts) (6.14) 
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is the covariance matrix between and y^s and 



(6.15) 



al^ (I^, ® C,3Cg) + E{(I^, ® AN|jAh,,AhJ,(I^, ® ANf^^} + ® I^,; 

(6.16) 



is the covariance matrix of y^s. The expectation in (16.161) is taken over all the random 
variables including A, Ah^^ f and lid,t of which the last two are functions of the random 
matrix H„. With the law of iterated expectations, i.e., E{X} = E{E{X|F}}, the second 
term of (I6.16P can be written as 

EhJE^,h,,{(I^. ® AN£;E{Ah,,AhJjH,,,Hj(I^, ® AN|;^|H^ (6.17) 

where a fact that the random matrix A is independent of Ah^^t is used. Since the term 
EjAh^^tAh^jHrf^i, H„} is not easy to tackle, we made an assumption that lid,t is Gaussian 
distributed under a given H^- In this case, AH^^ = — H^^ ^ is also Gaussian distributed 
and so does its vector form Ah^^ f. We know that AH^^ t is uncorrelated to H^ ^ refering to 
the orthogonality principle, therefore Ahd^t and H^^i are independent due to our imposed 
Gaussian assumption. The equation in fl6.17p 48^fhen given by 



EfijE^g^ J(I^, ® AN£;E{Ah,,AhJ|Hj(I^, ® AN|)^|HJ} 
Substituting (16.101) and the fact that Ng = Ijy^^^L into (I6.18p . we obtain 



(6.18) 



{N, - NL)al 



-E< 



-1 



®I7V, (6.19) 



The Hermitian term H^jH^ can be factorized into 



H^^f = UAU 



H 



(6.20) 



where U G C^^^^^ is the matrix whose columns are consisting of eigenvectors of H^H^ 
and A = diag(Ai, . . . , Aat^) is the diagonal matrix with diagonal elements being nonzero and 
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unordered eigenvalues of HuH^. Since the elements of both the uplink channel H„ and the 
noise matrix W2 are i.i.d. Gaussian distributed and the reverse training is assumed to 
be semi-unitary, each entry of Hu is i.i.d. Gaussian distributed. Referring to [10], we know 
that H^H^ has a Wishart distribution with A^^ degrees of freedom and its mean is given by 



E{H„Hf } = Nt 



where 



a 



2 A 



(6.21) 



is the variance of each i.i.d. random variable of H^j. Both the density function of HuH^ 
[TU] and the Jacobian of the eigenvalue value decomposition of H„H^ [TT] can be divided 
into the product of functions of A and U, thus we conclude that A and U are independent. 
With the independency and applying the law of iterated expectations, the equation in fl6.19p 



becomes 

{N, - NL)al 



ail 



hM. - j U . E.{ (/3A- + i.o-}u 



H 



-E> 



Nt 



(6.22) 
(6.23) 



./3/Ai + l 

where the equality holds since the eigenvalues of the Wishart distributed matrix H^H^ 
have identical distributions as any one of the unordered eigenvalues [12]. Replacing fl6.23p 
in f l6.16|) . we have an approximation of the covariance matrix of y^s as 



(AT, - NL)al a 



4 



a 



/3/Ai + 1 



+ < 



(6.24) 
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The normalized mean squared error (NMSE) of Hd can be computed as 

Tr(E{Ah,Ah^}) 

"^hj^VLS "^VLiVL-i ^ hay Li ) 



NMSE/ 



NtNL 



(6.25) 
(6.26) 



NtNL 

Substituting f l6.14p and fl6.24p into (16.261) . we have an approximation for the NMSE of the 
LR as 

/ 

NMSEl ^ 



1 



N 



\ 



(iV, - NL)al a 



'Ha 



a 



1 



J 

(6.27) 



For A^t ^ 1, the distribution of the eigenvalues of H„H^ is asymptotically approximated to a 
Gaussian distribution [13j, that is Ai ^ Af{Nta'^, Ntcr^) where is given in (I6.2ip . However, 
the expectation term in fl6.27p is intractable, we instead apply the Jensen's inequality and 
take its lower bound as an approximation. Hence, we have 



NMSE, 



1 ^3 
h 




a 



Nt 



\ 



(AT, - NL)<yl a 



a 



Nta 



7 



(6.28) 



On the other hand, the NMSE performance of the UR is analyzed as follows. The received 
signal of UR (15. 8p can be vectorized as 



yu-i = yec{Yu3) = {Inu ® Cts)g + (I^^ ® AN| )g + vg 



(6.29) 



where Cj3 = ^J~^Ct3, g = vec(G) and V3 = vec(V3). The covariance matrix between g and 
yu3 and the covariance matrix of yj/a are respectively given by 



^g,yu3 — ^gC^Nu ® C^) 



(6.30) 
(6.31) 
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Hence, the NMSE of the UR is given by 
^^^^^ = NJ^u 



NtNu 

al ' N,aUN,-NL)al + al 



Z2+4 .2r.r N,. , ,J (6-32) 



6.3 Optimal Power Allocation between Training and 
AN Signals 

With fl6.28p and (16 .32 p . we can jointly design the power values of {Sq, £1, £2, £3, cTq} 
considering the following power allocation problem 

min NMSEl ■■ '.■<ri^^"ir%: ■ (6.33a) 

s.t. NMSE[7 > 7 ifef§Ttii|J^^!i ^ (6.33b) 

£0 + £1 + £2 + £3 + (#;^|^aJiVj < i^^^(A/^ + iV^ + + iV^) (6.33c) 

£o + £3 + {Nt - NL)alNt < "Pt{Nt + Nt) (6.33d) 

£i+£2<PL{Nt + NL). (6.33e) 

Here, we aim to minimize the NMSE of LR subject to the constraint that the NMSE of UR 
remains above a preset lower limit 7. We also consider the average power constraint Pave 
and two individual power constraints Pf and P^ at the transmitter and LR, respectively. 
However, the problem is not easily solvable. To obtain an efficient solution, we resort to 
the monomial approximation and the condensation method often adopted in the field of 
geometric programming (GP) [Hj. Details are given in the Appendix 19. 2[ 
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Chapter 7 

Numerical Results and Discussions 



In this chapter, we present numerical resuhs to demonstrate the effectiveness of the proposed 
DCE schemes. We consider the MIMO wireless system as described in Chapter [2] with 
Nt = 4, Nl = 2 and Nu = 2. The elements of the channel matrices H and G are i.i.d. 
complex Gaussian distributed with zero mean and unit variance {ajj = cr^ = 1). Each entry 
of additive white noise matrices W, W and V is also i.i.d. complex Gaussian distributed 
with zero mean and unit variance, i.e., a"^ = a^^ = = 1. Moreover, the training lengths 
are set to be the antenna number of the terminal which transmits that training signal, i.e., 
tr = Nl = 2 and Tp = Nt = 4 for the reciprocal casqj and tq = = Nt = 4 and T2 = Nl = 2 
for the non-reciprocal case. Note that the overall training time is larger than the sum of 
all training length due to the processing time at the transmitter. Besides, the individual 
power constraints of the transmitter and the LR are respectively assigned as Pt = 30 dB and 
Pl = 20 dB. We incorporate an NMSE lower bound for comparison. The lower bound for 
reciprocal and non-reciprocal case are respectively given by 

nmselb ^(J^^'-'-jP'^'-P-'f^^my p.l) 



^In pure channel estimation, it is preferred to keep the training length minimum in uncorrelated channel 
and white noise [S] . We show in Figure 17.11 that the training length is better to choose as smallest length 
i.e., the number of transmit antenna in the secrecy channel estimation scheme. 
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Figure 7.1: NMSE performance versus forward training Tp for Pt = 30 dB and Pl = 20 dB. 



and 



NMSElb 



nonrec \ ^1 



1_ min{2A iVt, Pa.e(3iVi + iVi)} 

2 + 



(7.2) 



which both stand for the minimum achievable NMSE at the LR when a\ = 0, i.e., no AN 
exists. 

Figure 17.11 shows the NMSE performance of LR versus the forward training length Tp 
under the constant energy constraints. In the reciprocal case, the average energy constraint 
is given by Pave{Nt + Np) and the individual energy constraints of the transmitter and LR 
are respectively given by PtNf and PpNp. We compare different average power constraints 
Pave = 15 dB, 20 dB and 25 dB and different lower limit values 7 = 0.1 and 0.03. We see 
from Fig. 17.11 that the NMSE value of LR is monotonically non- decreasing with respect to 
the training length Tp. It shows that in secrecy channel estimation it is better to keep the 
training length as small as possible. This is due to the fact that as the forward training 
length increases, it takes more AN energy to satisfy the lower limit constraint on the UR 
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Figure 7.2: Power allocation among reverse and forward training powers Sr/t^i, £f/'Tf and 
AN power (A^"* - A^l)(t2. 

thus the budget for the training energises |s sacrificed. Moreover, we see that the lines of 
7 = 0.1 are more steeper than those of ■7 i!^-;Qr.03;..-The trade-off between the AN energy and 
training energies is more explicit as the lower limit is severer. 

Figure 17.21 shows the optimal allocation of the reciprocal case among the reverse and 
forward training powers £r/tr, Ef/tf and the AN power {Nt — NL)a1 versus average power 
constraint Pave- We compare two different lower limit values 7 = 0.1 and 7 = 0.03. We see 
from Fig. 17.21 that it is desirable to allocate more power to the AN and less power to the 
forward training as 7 increase from 0.03 to 0.1. This is due to the fact that the forward 
training signal benefits the LR and the UR equally while the AN primarily degrades the 
UR's estimation performance. In addition, we see that the reverse training power increases 
with 7, since the reverse training power mainly determines the subspace into which the AN 
is transmitted, which helps to reduce the interference caused by the AN on the LR. Note 
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Figure 7.3: Power allocation among the tr'^fi^_gi"^^ers £q/ Nt, £1/ Nt £2/Nl £^/Nt and AN 
power (At - A^)cr^ ff'^\SM^'^^^''^^>-- " 

that when Pave = 10 dB and 7 = 0.03, .this,. is the case with 7 out of the interesting interval 
f l4.10p . the reverse training power and AN p'^j^r-^'oth equal to which can not be showed 
in the log- value. 

On the other hand. Figure 17.31 shows the power allocation of the non-reciprocal case 
among the round-trip training powers £q/ Nt and £i/Nt, reverse and forward training powers 
£2/Nl and £3/Nl and the AN power {Nt — NL)al versus average power constraint Pave- We 
have similar observation about the allocation between the forward training power and AN 
power to that of the reciprocal case. We see from Fig. 17.31 that the round-trip and reverse 
training powers all increase with respect to 7, since these powers play the role to design 
the placement of AN for minimizing the interference on the LR. As 7 increases, so does the 
AN power, there needs more round-trip and reverse training powers to decrease the damage 
cause by the AN on the LR. 
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Figure 7.4(a) and Figure 7.4(b) show the NMSE performance of the LR and UR versus 
average power constraint Pave respectively for the reciprocal and the non-reciprocal channel. 
We compare two different lower limit values 7 = 0.1 and 7 = 0.03 in both figures. From 
Fig. 17. 4^ we observed that the NMSE of the UR meets the lower limit in both reciprocal and 
non-reciprocal case. Furthermore, the proposed DCE scheme constrains the UR's NMSE 



well above 7. In addition, from Fig. 7.4(b) , we see that the approximation of LR's NMSE 
fl6.28p is quite close to the Monte-Carlo simulation result of LR's NMSE. 

In Figure 17. 5[ we show the symbol error rate (SER) at LR and UR versus the average 
power constraint Pave in the data transmission phase. We consider the scenario where the 
transmitter sends a 4 x 4 complex orthogonal STBC (OSTBC) with Nt = 4. The code 
length is equal to four and each code block contains three QAM source symbols The 
data transmission power is set to Pave- Both LR and the UR will exploit their channel 
estimates obtained by the proposed DCE -tGitfe^d^' the received symbols. In this Monte- 
Carlo simulation, the SER is obtained by averaging over 50000 channel realization and 



OSTBCs. In particular. Figure 7.5(a) presents the associated average SERs for 64-QAM 
OSTBC in the reciprocal case. We see that the SER of the LR will gradually improve while 
the SER of the UR remains larger than 0.1 due to the poor channel estimation performance 



at the UR. Figure 7.5(b) shows the associated average SERs for 64-QAM OSTBC in the 
non-reciprocal case. We have similar observation in this case. Both figures illustrate that, 
with the proposed two-way training DCE scheme, the discrimination of the data detection 
performances between LR and UR can be effectively achieved. It is worthwhile to mention 
that the feedback-and-retraining DCE scheme proposed in [3j assumes a perfect feedback 
channel with no power consumption and, thus, it is difficult to have a fair performance 
comparison between the proposed scheme and that in [3j. 
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Figure 7.4: NMSE performance of the proposed DCE scheme for the reciprocal and non- 
reciprocal case. 
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Figure 7.5: SER performance of the LR and UR in an OSTBC system with the channel 
estimates obtained by the proposed DCE scheme. 
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Chapter 8 



Conclusion 



In this thesis, we proposed a new DCE scheme based on the two-way training methodology, 
where both the transmitter and LR are allowed to emit training signals. In particular, 
training signals sent by LR are used to help the transmitter obtain an accurate estimate 
of the transmitter-to-LR channel. The proposed two-way DCE scheme utilizes two phases 
of training in reciprocal channels and three phases of training in non-reciprocal channels. 
The proposed training design drastical,ly decreases the overall training overhead compared 



to the original DCE scheme proposed i]ST^:^'3?]M&^^aihing and AN powers were optimized by 
minimizing the NMSE of LR subject to a preset lower limit on the NMSE of UR, an average 
total power constraint, and individual power constraints over all transmitters. For the case 
with reciprocal channels, the optimal power allocation problem was reformulated into a one- 
variable optimization problem which can be easily solved by simple line search. For the 
case with non-reciprocal channels, we derived an approximation of LR's NMSE and utilized 
monomial approximation and condensation method to obtain an approximate solution for 
the power allocation problem. Numerical results were provided to verify the efficiency of the 
proposed two-way DCE schemes. 
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Chapter 9 
Appendix 



9.1 Proof of Proposition I 

For notational simplicity, let us define a = {Nt — NLja^. In the following, we solve the 
optimization problem in two steps: (i) find the optimal values of £p and a for any given £ji\ 
and (ii) find the optimal value of £r. .. //^f^^-^^. . 

Step i : i. ItetSiISi 

Suppose a feasible £r is given, the 6pCTQiAl!llb3flaes of £f and a can be found as functions 
of £ii from the below optimization problem 

{NLaj + aj,£R)£F 

max 2 (9.1a) 

s,,a>o NLal + al-£R + NLay^-a 

s-t- 2' ' 2 < 7, (9.1b) 

£f ^ ■ rp <Pave{TR^'rp) - £r, (9.1c) 

£f + a ■ tf < PtTF- (9. Id) 

Note that from (14 .Sp a feasible £r must satisfy £r < PlTr. In the following, we consider two 
different ranges of £r. 

Case 1 {PgveirR + tf) ~ J < £r < PlTr)- If Pavei^R + Tf) - ^ < PlTr holds. Since the 
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objective function in (19. ip is monotonically increasing with respect to Ep but decreasing 
with respect to a, by fl4.12p and the condition of Er > Pavei.'^R + 'Tf) — 1, we get £p{£r) = 
Pave{TR + Tp) — £r, a* = and hence the value of (19. lap becomes 

£;(^r) = PaUrR + Tp) - £r, (9.2) 

which is less than 7. 

Case 2 (Sr < minjPiTR, PaveijR + tf) — 7}): It can be observed that if the constraint (]9.1bp 
is inactive we can always decrease a until activating the constraint to obtain a larger ob- 
jective value. If the condition fl9.1bp is still inactive even when a = 0, we can instead lift 



8p to achieve a larger objective value while still satisfying (14. 9p . (I4.12p and the condition 
£r < ininlPiTR, Pavei^R + Tp) — 7}. We conclude that constraint f l9.1bp must be active at 
the optimum. Hence we have 



£^{£R)^/fft'^^ct{ER) + l]. (9.3) 



<:.. i 



By substituting (19. 3p into (19.11) . the problem becomes.; 

max + (9.4a) 

S.t. {Tp + ^^a + £R<Pa.e{TR + Tp)-^ (9.4b) 

7 {^01 + l^+Tp-a< PtTp. (9.4c) 
The range of Sr in this case is further divided into the following two subranges. 



(a) When Sr < ji and Sr < min{PLrjj, Pavei^R + Tp) - 7}, where ^ = Npi ^-S^ , 
the objective function in (I9.4ap is a monotonically decreasing function with respect to 
a. Therefore, the optimal value of a*{SR) is and the corresponding optimal objective 
value is equal to 7. 
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(b) When n < £r < min{PLr/j, PaveijR + tf) — 7}, the objective function in f l9.4ap is 
monotonically non- decreasing with respect to a. For ii < £r< Pavei'TR + Tp) — PtTp, 
constraint f l9.4cp must be active at the optimum with 

Tp + a'c^ al 



Reversely, considering £r > max{/i, PaveiTR + Tp) — PtTp}, the constraint f l9.4bp must 
be active at the optimum with 

N Pave{TR + Tp) - ^ - Sr 

a [£r) = ■ — (9.6) 

Moreover, for £r > fi, the optimal objective value of (I9.4ap can be shown to be 

^hc*/a'^ + 1 



which is no less than 7. 
Step ii : 



■ 7 > 7- (9.7) 



We now solve for the optimal value of Sr. From the analysis in Step i, a feasible Sr 
satisfying Sr < mm{PLTR, Pavei^R + Tp) — 7} leads to greater objective value than that of 
Pave{TR + Tp)—'^ < Sr, < PlTr, thus the Optimal value of Sr must lie in the former condition. 
For the first case that /i > mm{PLTR, Pave{TR + Tp) —7}, we can infer Sr< ji for all feasible 
Sr satisfying Sr < mm{PLTR, PaveijR + Tp) — 7} so that a* = and thus Sp = 7. Then we 
get S^ = for no need of AN. For the other case of /i < mm{PLTR, PaveijR + Tp) — 7}, from 
(19. 2p and (19.71) we can see that the corresponding objective value for max{0, ^i] < Sr < 

TSim.{PLTR, Pave{TR + Tp)-^} is UO leSS than that for Sr < H. If H < Sr < Pave{TR + Tp)-PtTp 

exists, the optimization problem (19. 4p becomes 

max iNp^l + ^'HSR)S^.^ 

S.t max{0, H} <Sr< PaveijR + Tp) - PtTp 
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where Sp and a* are given by (19. 3p and (19. 5p which do not depend on £r in this condition. 
It can be observed that the objective function (19.81) is monotonically non- decreasing with 
respect to a; thus the optimal value is achieved when S'^ = Pavei^R + Tp) — PtTp. However, 
the corresponding optimal objective value of (19. 8p is the same as the objective value of (I4.13P 
in this case. Hence, we can have the value of S'^ lie in the interval max{0, fi, Pavei'TR + 'Tp) ~ 
PtTp} < ^^_R < HiinlPiTij, PaveiTR + Tp) —7} by solviug the optimization problem (I4.13P and 
the corresponding a{Sfi) and £p{£ji) are given by (19. 6 p and (19. 3p . respectively. 

9.2 Monomial Approximation and Condensation Method 
for the Problem in (16.331) 



Here, we show how to obtain an efficient solution for the problem in (I6.33P using monomial 
approximation and condensation method. The problem in (I6.33p can be stated as follows: 

/II /l(go,C.^g2,g3)V ' . . 

mm ^- + ^-TT? 1 c c \ (9-9a) 



subject to TT ^777 ^r-^ t < (9.9b) 

8^ + 81 + 82 + 83 + {Nt - NL)alNt < PaveiNt + Nt + NL + Nt) (9.9c) 

80 + 83+ (Nt - NL)alNt < Pt{Nt + Nt) (9.9d) 

81 + 82 < P^Nt + Nl) (9.9e) 



where 
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and 



3 



^ I I ^lf£^^^2\ 2 , ^^1*^ /^^|/o , 2\ 2^2 , ^2 , 

+ I I \- Cr,„ I tt H ^ I \- Cr,„ I tt — + — + 



(9.10) 



and by defining the variables to = -7^ (^ti ti = ^ h = ^3 = j^, and ^4 



Nt ■ ■ V V AT,. 0-1,^ 

By introducing the auxihary variable 

^ ^ fi{So,a^,£2,£3) 
/2(£^o,a^^^2,^^3) 

{Nt — NL)alaQ + a^, the problem can be reformulated as 

min t"^ (9.11a) 
subject to t < ^f°'^^'^^'^^| ,..wv,. (9.11b) 

citatr^ < 1 ^^^-^Mf (9.11e) 
C2 -^to + NtNMi + NLt2 + A^ttg + -^U < 1 (9.1 If) 

cs ( 4^to + Nth + ) < 1 (9.11g) 
Ci{NtNLUi + NLt2)<l (9.1 Ih) 



where 



fl{to-,tl,t2,t3) —Nitotits H ir-^tQtit2H + ^2^3 H ^^3; 

^1 



/t4 










kL( 


V4 







/2(t0, tl, t2,h) -[^-^] (Th, ( -f^^O^l + * 2^" *1^2 + -^t2 + 2 ^2 



4 



^S, (^w ^Hu^w 

2 ^2 



+ iVLto^l + f^totlt2 +t2 + 
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and 



Cl = I 2" J 5 C2 = I Pavei^Nt + Nl) H 3 ^ 2~ 

C3 = ( 2PM + ^ + C4 = (P.(iV, + iV,))-^ . 

V '^H^ ^G / 



Note that in f lQ.llbp the equahty was replaced by the inequahty since one can inspect that 



the inequahty must be active when the optimal objective value is achieved. To make sure 
and a\ are no less than zero, we attach two posynomial constraints (19.1 Id) and (IQ.lldp . In 
addition, we have reformulated the constraints in f l9.11etl9^11hp into posynomial inequalities. 



which are standard inequality constraints for GP. However, the inequality constraint in 
flQ.llbl) is not a standard GP inequality. It can only be expressed as a ratio of posynomials 
as given below: 

i^tQtit^t H -^tit2t4t + + -jT^t^t) + Nit^tit H -jf^tQtit2t + t2t + -jT^t 



< 1 



w Hu 



3 



(9.12) 

In order to simplify the problem into a. standard GP form, we apply the monomial approxi- 
mation fT^ to transform this into a posynomiaf constraint. In particular, if a set of feasible 
points {t, ti, ^2, ^3, ti} of problem (19. lip is given, the inequality in (I9.12p can be replaced by 
the posynomial constraint given below: 



(^totit4t H -J^tit2t4t + -\-t2t/\t + s 2 ^4i^) + Nit(jtit H -2^^0^1^2^ + ^2^ + 'ZT^'^ 



gitMMM (I)" (l)'^ (l)'^ (l)'^ 



< 1 



(9.13) 



where 



-'^d / 4- 4- I J^u J. J. I J. I 



^(t, tl, ^2, ts) =^^^ ^^0^1 + f^tit2 + ^t2 + 



+ Nitotits H 7r-^totit2t3 + ^2^3 H 5— '^3, 

^5) 
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and 

Hence, for a given set of feasible points {t, ti, ^2, ^3, ^4}, the problem in (19. lip can be approx- 
imate by the following problem 

min (9.14a) 

eo, 61,62, 63, (t2>0 



subject to 



\ '^u' '^iB '^■ui '^Hu ™ 



^(Mi,t2,t3)(|) (f^) (1) (1) 

cr^t-^ < 1 (9.14c) 

< 1 (9.14d) 

cihtl^ < 1 (9.14e) 

C2 -^to + NtNMi + NLt2 + Nth + -4U < 1 (9.14f) 

ca + Nth + < 1 9.14g 

c^{NtNLhh + NLt2) < 1 (9.14h) 

The problem then becomes a standard GP problem and can be efficiently solved by a few 
simple computer softwares such as CVX [15]. The condensation method then proposes to 
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repeat this process iteratively by replacing the set of feasible points {t, ti, ^2, ^3, t^} in each 
iteration with the optimal solution of f l9.14p obtained in the previous iteration. This process 
continues until no further improvements can be obtained. 
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